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ABSTRACT 

There is now compelling evidence of a link between long-duration gamma-ray 
bursts (GRBs) and Type Ib/c supernovae (SNe). These core-collapse explosions 
are conjectured to radiate an anisotropic, beamed component associated with a 
decelerating, relativistic outflow and an unbeamed, isotropic component associ- 
ated with the slowly expanding stellar debris. The anisotropic emission remains 
at a very low level until the Doppler cone of the beam intersects the observer's 
line of sight, making off-axis GRB jets directly detectable only at long wave- 
lengths and late times (> 1 yr). Circumstellar material, however, will Compton 
scatter the prompt gamma-ray and afterglow radiation flux and give rise to a 
reflection echo. We show that the Compton echo of a misaligned GRB carries an 
X-ray luminosity that may exceed by many orders of magnitude that produced 
by the underlying subrelativistic SN during the first few weeks. Bright scattering 
echoes may therefore provide a means for detecting a population of misaligned 
GRBs associated with nearby Type Ib/c SNe and yield crucial information on 
the environment surrounding a massive star at the time of its death. The ques- 
tion of whether the interpretation of GRB980425 as an ordinary GRB observed 
off-axis is consistent with the lack of an X-ray echo is addressed, along with the 
constraints derived on the possible existence of misaligned GRB jets in SN1993J, 
SN1994I, SN1999em, and SN2002ap. 

Subject headings: gamma rays: bursts - gamma rays:theory - supernovae: gen- 
eral - X-rays: general- stars: Wolf-Rayet 
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1. Introduction 

The first hint of a connection between core-collapse supernovae (SNe) and gamma-ray 
bursts (GRBs) came in April of 1998, when SN1998bw was detected in the BeppoSAX error 
box of GRB980425 (Galama et al. 1998). While unremarkable in its apparent properties 
(peak flux, duration, burst profile), GRB980425 (isotropic) energy release of only 8 x 10 47 
erg was some four orders of magnitude smaller than that of typical GRBs (e.g., Bloom et 
al. 1998). The observational basis for a connection between GRBs and Type Ib/c SNe 
has been greatly strengthened by the recent discovery (Hjorth et al. 2003; Stanek et al. 
2003) of a SN1998bw-like spectrum lurking beneath the optical afterglow of GRB030329. 
This SN has been designated SN2003dh. Due to its extreme brightness and slow decay, 
spectroscopic observations of GRB030329 have been extensive. The early spectra consisted 
of a power-law decay continuum (F v oc z/ -0 9 ) typical of GRB afterglows with narrow emission 
features identifiable as Ha, [OIII], H/5 and [Oil] at z — 0.1687 (Greiner et al. 2003), making 
GRB030329 the second nearest burst overall (GRB980425 being possibly the nearest at z — 
0.0085) and the classical burst with the lowest known redshift. At this distance, GRB030329 
is typical in its gamma-ray budget - its total isotropic energy release is 9 x 10 51 erg (Hjorth 
et al. 2003). 

Independent of their connection with a GRB, SN1998bw and SN2003dh belong to the 
unusual class of Type Ic SNe. The lack of hydrogen lines in both spectra is consistent with 
model expectations that the stellar progenitor lost its hydrogen envelope to become a Wolf- 
Rayet (WR) star before exploding (Woosley 1993). The broad lines are also suggestive of 
an asymmetric explosion viewed along the axis of most rapid expansion (e.g., MacFadyen & 
Woosley 1999). The radiated energy may therefore be considerably smaller than the total 
explosive energy - in the range 1 — 2 x 10 52 erg - inferred assuming spherical symmetry 
(Woosley et al. 1999). Moreover, the larger ratio of afterglow luminosity to SN luminosity 
in GRB030329 is consistent with a burst that has been observed nearly on axis, whereas 
the low luminosity of GRB980425 and its afterglow hint to a substantially off-axis view 
(Nakamura 1999; Eichler & Levinson 1999; Woosley et al. 1999; Granot et al. 2002). 
In these models, together with the afterglow radiation restricted to a narrow relativistic 
beam, there is also quasi-spherical emission from slower moving SN ejecta running into the 
circumburst medium (Ramirez- Ruiz et al. 2002; Zhang et al. 2004). The flux associated 
with the beamed component declines rapidly with observer direction, thus making the task 
of detecting misaligned GRB jets particularly challenging (Granot et al. 2002). 

Even if GRBs were caused by a standard set of objects, their appearance would depend 
drastically on orientation relative to the line of sight. The gamma-rays we receive come 
from material the motion of which is directed preferentially along the rotation axis of the 
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stellar progenitor, and therefore provide no information about the ejecta in other directions. 
The propagation of gamma-rays through an external medium would, however, randomize 
the initial photon pulse energy and, because of light-travel time effects, produce a luminous 
reflection echo (Dermer et al 1991; Madau et al. 2000; Ramirez-Ruiz et al. 2001; Sazonov & 
Sunyaev 2003). We show below that the X-ray luminosity of such an echo exceeds during the 
first few weeks that produced by the underlying subrelativistic SN. The detection of energetic 
scattered light, with luminosities high enough to be detected by Chandra and XMM-Newton 
observations (and the Swift experiment now under construction), may therefore provide a 
means for uncovering a population of misaligned GRB jets associated with nearby Type Ib/c 
SNe. Conversely, one may use the lack of evidence for a Compton echo to set constraints on 
the the density profile of the circumburst material that surrounds a massive star at the time 
of its death. 



The deceleration of a relativistic jet propagating through a stellar wind starts when 
about half the initial energy has been transferred to the shocked matter, i.e. when an 
amount Tq 1 of its own rest mass has been swept up. The relativistic expansion is then 
gradually slowed down, and the blast wave evolves in a self-similar manner with a power- 
law lightcurve (Blandford & McKee 1976). This phase ends when a total mass E n /c 2 (we 
denote with Eq the burst energy per unit solid angle) shares the energy and the Lorentz 
factor drops to 1. As, in the unshocked stellar wind, the mass within radius r is Mr/v w , the 
blast wave deceleration radius is rd cc = Eqv w / \Mc 2 Tl) k (2 x 10 15 cm) Eq^Vw^MZ^Tq^, 
where M_5 = (M/10 -5 M Q yr _1 ) and v w ^ = (w U) /10 3 kms _1 ). At observer times of more 
than a day, the blast wave would be decelerated to a moderate Lorentz factor, irrespective 
of the initial value (Rhoads 1999). The beaming and aberration effects are thereafter less 
extreme, so we observe afterglow emission not just from material moving almost directly 
towards us, but from a wider range of angles. A temporal break occurs when the angle 
9 ~ T^ 1 of the Doppler beaming cone equals the jet opening angle 6 - For a blast-wave 
decelerating according to the approximation T(r) = T for r < r dec , and T(r) = r (r/r dcc )~ s 
for r > rdec, the time tb at which the temporal break is observed after the explosion is given 
by t b (1 + z)[— r cos# + Jq dr'fi _1 (r)]/c. Here c(3(r) = c[l — T(r) -2 ] 1 / 2 is the speed of the 
blast wave. This gives 



2. X-ray Emission from Misaligned GRB Jets 



2.1. Off-axis Afterglow Emission 
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The lightcurve seen by an observer located within the initial jet aperture 9 < 9 is very 
similar to that for an on-axis observer (e.g., Dermer et al. 2000) but differs markedly in 
intensity once the observer is outside the beaming cone. The emission remains at a very 
low level until the Doppler cone of the beam intersects the observer's line of sight. This 
can be seen by comparing the 9q and 29q curves in Figure 1. Following the arguments given 
to derive equation (1), we find that an off-axis observer {9 > 9 ) starts to detect emission 
at a level comparable to an on-axis observer at times given by equation (1), but with 9 
replaced by 9 — 9 . As shown in Figure 1, it appears that the X-ray emission from highly 
misaligned GRB jets (i.e. 9 > 39 ) is not significantly larger than that emitted by the SN 
remnant itself. Since the jet opening angle of a GRB is typically assumed to be of the 
order of 9q < 10° (Frail et al. 2001; which represents 0.76% of the full sky) it is fair to 
conclude that detection of short wavelength transients associated with core-collapse SNe will 
not unambiguously demonstrate the presence of a GRB jet. The ratios of the on-axis to 
off-axis fluxes are however much smaller at radio frequencies where the emission persists at 
comparable levels for about the same length of time (Paczyhski 2001). 



The beaming angle for the gamma-ray emission could be far smaller than that of the 
late-time afterglow, and is much harder to constrain directly. If the gamma rays were much 
more narrowly beamed than the optical afterglow there should be many orphan afterglows. 
The sudden brilliance of a GRB propagating in a preburst stellar wind will be reflected by 
the circumstellar gas. If E = J E&dedtt is the total energy emitted by the burst, where E € q 
is the energy emitted per unit energy e and unit solid angle Q, then the equivalent isotropic 
luminosity of the Compton echo inferred by a distant observer is 



where r is the distance from the site of the burst, da/dQ is the differential Klein-Nishina cross 
section for unpolarized incident radiation, e' = e[l + e(l — cosy?) /m e c 2 ] -1 is the energy of the 
scattered photon, and n w is the local electron density (Madau et al. 2000). The equal-arrival 
time scattering material lies on the paraboloid determined by r = ct(l — cos(/?)~ 1 /(l + z), 
where cp is the angle between the line of sight and the direction of the reflecting gas as seen by 
the burst (e.g., Blandford & Rees 1972), t is measured since the burst is first detected, and 
{1+z) is the cosmo logical time dilation factor for a source at redshift z. At energies below 100 
keV, where the echo mirrors the spectral energy distribution of the prompt pulse, echoes will 
typically be harder than the X-ray afterglows observed at around 1 day. Above 200 keV, and 
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Reflection Echoes 
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independently of the incident gamma-ray spectrum, reflection echoes will show a high-energy 
cutoff because of Compton downscattering (Madau et al. 2000). One distinctive feature of 
the echo signal is that the scattered radiation is likely to be polarized. Even when the incident 
beam is unpolarized, the expected fractional polarization is II = (1 — cos 2 (/?)/(l + cos 2 y?). 
Therefore, in the case of observing the jet precisely at ip = 90°, the scattered radiation is 
100% polarized. On the other hand, if the prompt radiation is 100% polarized (e.g., Coburn 
& Boggs 2003), no net intensity will be seen at (p — 90°, since the electron's motion is 
confined to a plane normal to the observer 2 . 

The reflected flash luminosity of a collimated GRB propagating through a constant 
velocity wind with n w (r) = Ar~ 2 is dominated by the receding jet as, at a given observer 
time, it originates closer to the GRB where the density is higher. In this case, 

<** = 3 -^§T L = ( 3 x 1Q41er S O ^)M_ 5 < 3 £ 5 i (1 + zf t~l, (3) 

(Madau et al. 2000), where ot is the Thomson cross section £51 = i?/10 51 erg is total energy 
of the collimated GRB pulse, and ip(<p) = (1 + cos 2 ip). The echo declines with time faster 
than most afterglows observed on-axis but exceeds them in intensity once the observer is 
outside the beaming cone (Fig. 1). For this reason, backscattered radiation could provide 
a means for detecting a population of misaligned GRBs in nearby core-collapse SNe and 
provide important information about the circumburst environment. 

Equations (2) and (3) assume that photons scatter only once and that absorption is negli- 
gible. As the reflecting zone propagates along the receding beam, the backscattered radiation 
on the way to the observer will pass with impact parameter b = r sin ip through regions of even 
higher density than the scattering layer. For a jet inclined at angle p> to the line of sight, the 
electron column density traversed by the backscattered radiation in the steady, spherically 
symmetric wind of a WR star is N w = A(n—ip>)/b « (6x 10 19 cm -2 ) ^{<p) M_ 5 v~^ 3 (1+z) , 
where ^(<p) = (1 + cos <^) (7r — ip) / simp. This amounts to an optical depth for Thomson scat- 
tering of r T ~ [10~ 2 , 2 x 10~ 4 ] t da y S (l + z) for ip = [1°, 45°]. Photoabsorption at 1 keV takes 
place with a cross section of ~ 2 x 10~ 22 cm 2 (for solar abundances). Multiple scattering 
and photoabsorption can therefore be neglected all times 3 but in the very early light curve. 



2 Thc differential cross sections for 100% polarized and unpolarized radiation, are in the ratio 2 cos 2 <p : 
1 + cos 2 (p (e.g., Rybicki & Lightman 1979) 

3 The optical depth to electron scattering could in principle be much higher if some lines of sight to 
the observer traversed the SN ejecta. This is unlikely, however, as it requires wsn/c > sin< f o/(l + cost/?), 
corresponding to wsn/c > 0.2 for jets with ip > 20°, where i>sn is the velocity of the SN ejecta. 
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3. A GRB Behind the Veil? 

Collimated GRB outflows directed away from an observer will produce extended X-ray 
transients at a flux level orders of magnitude lower than transients observed along the jet axis. 
The off-axis jet interpretation of GRB980425 requires the viewing angle to have been > 39 
with 9 the half angle of the most energetic part of the jet, otherwise the X-ray afterglow 
would have contaminated the SN light curve unacceptably (Fig. 1; similar requirements have 
been derived for the optical light curve by Granot et al. 2002). The off-axis jet interpretation 
of GRB980425 typically assumes that our line of sight is a few degrees from the sharp-edge 
conical jet (see dashed lines in Fig. 1). This interpretation would be inconsistent with radio 
observations since in this case the jet deceleration would bring our line of sight into the jet's 
radiation beaming cone, leading to strong emission on timescales < 1 year (Waxman 2004a). 
Our line of sight must make a large angle with the jet axis (> 49) in order to avoid observing 
strong radio emission. The observed gamma-ray flux of GRB980425 can be explained in 
this case by assuming that the jet is not sharp-edged, but rather has wings of lower Lorentz 
factors that extend to large 9 (Woosley et al. 1999) and produce the observed faint burst. 
Interestingly, this material would also increase the afterglow luminosity at early times when 
almost nothing is seen of the inner jet. At early times, the X-ray emission of the underlying 
subrelativistic SN is probably small (Waxman 2004b). The X-ray emission during the first 
few weeks should be therefore dominated by the rapidly fading echo (Fig. 1). The fact that 
the echo decay is both so steep and presents little spectral evolution will unambiguously 
distinguish it from any relativistic blast wave in which the electrons emit a constant fraction 
of the energy gained in the shock (Cohen et al. 1998). The lack of a detectable Compton 
echo therefore places severe constraints on the brightness of a possible GRB jet associated 
with SN1998bw, and it is to this problem that we now turn our attention. 

During the first day after the explosion no significant flux was observed by BeppoSax in 
the direction of GRB980425 to a limit of about 10 41 erg s _1 . An off-axis jet with standard 10 51 
erg energy expanding into a wind with M xs 10~ 5 M Q yr _1 and v w fa 1000 km s _1 can therefore 
be ruled out. A misaligned jet interpretation for GRB980425 may still be consistent with 
the observations if it is assumed that the wind of the progenitor was atypically-low density, 
with M_5/f Wj 3 < 0.1 (Fig. 2). The lower value of M^/v w ^ is consistent with the observed 
radio emission from the shock driven into the progenitor wind by the SN ejecta (Waxman & 
Loeb 1999; Li & Chevalier 1999), although, as usually in the case of radio SNe, the data are 
unable to brake the degeneracy between M/v w and e# (i.e. the fraction of the SN post-shock 
thermal energy carried by magnetic field). A model with es <C 1 and M_ 5 /v w ^ > 0.1 would 
need e# < 10 -4 in order to reconcile radio observations with the prediction of an off-axis 
jet (Waxman 2004a; Soderberg et al. 2004). If, on the other hand, one adopts a model 
with €b near equipartition, then M_ 5 /v W)3 < 0.1 is required. Our analysis demonstrates that 



-7- 



the lack of detection of a GRB echo places robust constraints - i.e. independent of e B and 
complementary to those derived via radio observations - on the assumed density profile. By 
the same token, similar constraints on the intrinsic energy of a misaligned GRB, compatible 
with the circumstellar density profile inferred from radio observations, can be placed for 
SN1993J, SN1999em, SN2002ap and SN1994I (Fig. 2). 

4. Summary 

GRB explosions are conjectured to radiate an anisotropic, beamed component from 
the decelerating, ultrarelativistic outflow and an unbeamed, isotropic component from the 
slowly expanding stellar debris. The flux associated with the former depends upon the 
observer direction and declines rapidly with observer time (Fig. 1). Hence the X-ray flux 
detected by an on-axis observer is much larger than the flux that would be detected from an 
identical source viewed off-axis. The relative flux levels of on-axis and misaligned jet sources 
are, however, much smaller in the radio/infrared than at optical and X-ray frequencies. 
Transients from misaligned GRB outflows could then be detected at late times with radio 
follow-up observations (e.g., Granot & Loeb 2003). Here we have considered an alternative 
possibility: the detection of scattered radiation from a narrowly collimated but misaligned 
pulse of gamma-rays propagating through a WR-type wind. While emission from the SN 
shock would dominate over the scattered radiation at late times because of the steeper 
lightcurve of a reflection echo, this will not be true at early times when the echo may be 3 
orders of magnitude brighter. Compton echoes in nearby SNe could be studied with Chandra 
and XMM-Newton observatories out to a distance of 0.1 Gpc. The planned Swift experiment 
may also soon provide some constraints for individual bursts at higher threshold energies. 

We have argued that the lack of a detectable echo in GRB980425 may be consistent with 
an off-axis jet interpretation only if the density of the circumburst wind is M_ 5 /t> W)3 < 0.1, 
lower than expected for a WR star. The mass-loss history of a WR star during its last 
few centuries could be quite complicated as the star enters advanced burning stages unlike 
those in any WR star observed so far. The lightcurve of the resulting echo will depend fairly 
strongly on the properties of the progenitor system, especially the mass-loss rate of the star. 
The deceleration of a presupernova wind by the pressure of the surrounding medium could 
create circumstellar shells, as would the interaction of fast and slow winds from massive 
stars. The complex density structure seen in SN1987A may be a hint of the immediate 
neighborhood of a GRB. Further data on X-ray emission from nearby core-collapse SNe will 
therefore offer important clues to the nature of the precursor star. 
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Fig. 1. — X-ray emission from misaligned GRB jets. Compilation of observed X-ray lightcurves for GRBs 
and core-collapse SNc (Kouveliotou et al. 2004) presented in the source frame using a cosmology with 
TO = 0.27, = 0.73, and Hq = 72 krns" 1 Mpc -1 . Dashed lines: Lightcurves calculated at various 
inclination angles 9 for a GRB with standard parameters: i?n,54 = 1, To = 100, p = 2.5, e e = 0.5, eg = 10~ 3 
and #o = 5° (where e# and e e are the fraction of the internal energy in the magnetic field and electrons, 
respectively, and p is the power- law index of the electron energy distribution). Solid lines: Compton echo of 
a GRB. Lower two curves assume L = L a tj ays with L a — 10 40 and 10 42 erg s" 1 (see eq. [3]). The primary 
burst is assumed to be a two-sided collimated pulse propagating at an angle ip ~ 1 rad and invisible to the 
off-axis observer. The upper curve (highest luminosity echo) is calculated for the dense environment expected 
at the end of the evolution of a 29 Mq WR star with solar metallicity (Garcia-Segura et al. 1996). The 
collimated GRB, whose total energy is 10 52 erg, is assumed to have a broken power-law spectrum, eE e oc e 
if e < 250 keV and e -0 - 25 otherwise. It is instructive to look at the contribution of the red-supergiant 
shell (of radius r s h) to the luminosity of the scattered echo. The time interval, 2r s h cos ip/c, between the two 
reflection spikes at times ~ 1000 — 2000 days can be used to determine <p. 



Fig. 2. — Constraints on the density of circumstellar material and burst energetics. The lack of an X-ray 
echo at about 1 day from the explosion poses severe constraints on the existence of a misaligned GRB jet 
in SN1998bw (shaded region). An off-axis jet interpretation may be consistent if the density of the stellar 
medium is lower than expected for the assumed WR progenitor. Constraints on the total kinetic energy 
derived from modeling of the radio emission of SN1998bw are shown as open circles (Waxman 2004a; Li & 
Chevalier 1999). Similar constraints on the possible existence of a misaligned GRB jet compatible with the 
circumstellar density profile inferred from radio observations (Weiler et al. 2002; black symbols) for various 
core-collapse supernovae are shown. Constraints derived by Panaitescu & Kumar (2001) from the modeling 
the GRB afterglow are also illustrated. 
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